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Abstract 

The exact fermion propagator in a classical time-dependent gauge field is derived by solving the 
equation of motion for the Dirac Green's functions. From the retarded propagator obtained in this 
way the momentum spectrum for the produced fermion pairs is calculated. Different approxima- 
tions and the exact solution for the propagator and the momentum spectrum are presented. 
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I. INTRODUCTION 



Particle production in classical bosonic fields has been a topic of continuing interest in 
quantum electrodynamics (QED) and quantum chromodynamics (QCD). It is relevant for 
the physics of the early universe of intense laser fields as well as of ultra- relativistic 
heavy- ion collisions and the quark-gluon plasma [3] (QGP). A lot of effort is made to study 
the QGP's production and equilibration JJ in nuclear collision experiments at the Rela- 
tivistic Heavy-Ion Collider (RHIC) at the Brookhaven National Laboratory (BNL) and the 
Large Hadron Collider (LHC) under construction at CERN. The existence of such a state 
of matter is predicted by lattice QCD calculations at high temperatures [Hf . 

At ultra-relativistic energies, the two nuclei are highly Lorentz contracted. When they 
pass through each other, a chromoelectric field is formed due to the exchange of soft gluons 
[6j]. This is a natural extension of the color flux-tube model or the string model which are 
widely applied to high-energy pp, e + e~, and pA collisions • Many other recent publications, 
e -g- |a liA llfl lU| are based on the hypothesis that the initial state in heavy-ion collisions is 



dominated by gluons which on account of the large occupation number can be treated as a 
classical background field. 

The larger the occupation number (c\c^ of the bosonic sector of a physical system, the 

better it can be described by a classical field. Here Ck and c k are the bosonic field annihilation 
and creation operators for particles of momentum k. Especially if the occupation number 
is much larger than one, the commutator of the creation and annihilation operator and 
hence quantum effects can be neglected: [c\, Cfc] = 1 (c\ck^. The field operators can be 
approximated by complex numbers, i.e. they are treated classically. 

For gluons in a heavy-ion collision at RHIC with = 130GeV the initial occupation 
number for gluons of average transverse momentum ~ lGeV in the center of the colli- 
sions is roughly equal to 1.5 2|- Although this number is not much larger than unity, the 



classical field as the expectation value of the gauge field still constitutes the main contribu- 
tion as compared to the fluctuations of the gauge field. The occupation number for lower 
transverse momenta is yet higher and the classical concept is an even better approximation. 
For larger transverse momenta the occupation number is smaller and thus the quantum 
fluctuations are more important there. As the above value for the occupation number at the 
average transverse momentum is larger than one, keeping only classical bosons is justified 
as a first approximation. One can investigate quantum fluctuations in a subsequent step. 

The high-occupation number bosonic fields are of the order A ~ g . Thus processes 
with multiple couplings to the classical field are not parametrically suppressed by powers of 
the coupling constant g. Without an additional scale, they have to be taken into account 
to all orders. Under the prerequisite of weak coupling, the most important quantum pro- 
cesses involve only terms of the classical action which are of second order in the quantum 
fields. These are the fermion and the antifermion fields as well as the field of the bosonic 
quantum fluctuations. The coefficient of the second order terms for a given field constitutes 
the inverse of the corresponding two-point Green's function. Inversion for selected boundary 
conditions yields the particle's propagator. The propagators contain all the information on 
two-(quantum-)particle reactions in the presence of classical fields to all orders in the cou- 
pling constant g. These reactions are scattering off the classical field or particle production 
by vacuum polarisation. 

In the following let us consider particle production. In quantum electrodynamics this 
means production of electron-positron pairs. Analogously, in quantum chromodynamics 
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quark- ant iquark pairs can be produced. However, due to the non-linearity of the field tensor 
pairs of gluonic quantum fluctuations are produced, too. In fields of the magnitude A ~ g^ 1 
the production of both kinds of pairs is equally parametrically favoured. This paper only 
deals with the production of fermions and antifermions. It is possible that in a given situation 
the bosonic sector is covered by the concept of a classical field sufficiently well. Corrections 
to the high- momentum sector could be perturbatively accessible fl3l |. 

No concept of a classical field exists for the fermionic sector. There the occupation 
number has always to be less than unity due to the Pauli principle. A purely perturbative 
treatment could only describe the high-momentum sector. The soft part would not be 
treated consistently. 

An alternative to the perturbative approach for particle production is Schwinger's 
constant-field method [lj] which is an exact one-loop non-perturbative approach. This 
method can also be understood as semiclassical tunneling across the mass gap |l5j |. How- 
ever, this scheme is based on the assumption of a slowly varying classical field. If the field 
changes too rapidly in space or time, the production of fermions is again not described prop- 
erly. A different concept is needed that is independent of energy or time scales, respectively. 
Such a concept is especially important, if the time scale for a process is to be determined 
by a self-consistent calculation. If an approach is applied for such an investigation, which 
relies already on an assumption about the time scale, the result is likely to be misleading. 
For example, if the decay time scale for a classical field is to be calculated based on particle 
production, perturbative concepts are likely to lead to times which are too short, while the 
Schwinger method tends to predict a development of the system that is too slow. 

For a concise treatment other methods are necessary. Exact results are desirable but hard 
to obtain. As mentioned before, neglecting bosonic quantum fluctuations the behaviour of 
the fermions is governed by their two-point functions in the classic background field. It 
can be obtained by solving the equation of motion for the Dirac Green's function G(x, y) 
exactly: 



There are other ways to derive the full propagator, for instance by resumming all terms 
of the perturbative series or by adding up a set of (at all times complete) wave-function 
solutions of the Dirac equation. 

In arbitrary fields a few general approximations are known. Neglecting the field in the 
equation of motion leads to the free Green's function G°(x — y). The standard perturbative 
series is a sum of terms containing powers of the background field A between free Green's 
functions. The asymptotic behaviour of the free Green's function determines that of the 
approximated full Green's function. Another approach which applies in an arbitrary field 
is the static approximation G s (x,y). It is obtained by neglecting the spatial part of the 
covariant derivative in the differential equation. The remaining ordinary differential equation 
can be solved by direct integration. Yet another approach can be found in [l6| . 

The following investigation concerns the case where the classical field depends arbitrarily 
on one rectilinear coordinate A = A(n-x). The equation of motion for the propagators shall 
be solved directly. If such solutions are investigated, care has to be taken: the result could 
be any Green's function which is not necessarily a propagator. If a propagator has been 
obtained, the imposed boundary conditions determine whether the result is the retarded or 
the Feynman propagator or one of their related singular functions. 




(1) 
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A solution in a field depending on one rectilinear coordinate can also be seen as an 
approximation for the case where the strongest dependence is on this rectilinear coordinate 
and the dependence on all the others is much weaker. 

It will be explored whether an approximation scheme can be found that is independent 
of assumptions on time and/or energy scales over a large range of parameters. The in- 
vestigated approximations are the Born, the weak-field, the strong-field, and the Abelian 
approximation. The weak-field approach is an expansion in powers of the gauge field A 
based on the free propagator and valid for A C w. w stands for the on-shell energy of the 
described particles. The strong- field approximation is justified for u <C A and consists of 
an expansion in powers of the on-shell energy. For the Abelian approach the commutators 
of the elements of the Clifford and the charge algebra are neglected. 

Chapter II includes the exact solution of the equation of motion for the Dirac Green's 
function and gives several approximations to the full solution. Section III contains the 
application of the previous findings to the problem of particle production and the comparison 
of the different schemes to the exact result. In the last chapter the contents of the paper 
are summarised. 

Throughout the paper the metric is: g^ u = diag(l, — 1, — 1, — 1), angular momenta are 
measured in units of fi, and velocities in fractions of the speed of light c. From hereon, the 
coupling constant is included in the classical field: gA% ld = A^ ew . 



II. DETERMINATION OF THE PROPAGATOR 

Let us consider homogeneous solutions Gn{x,y) of Dirac's equation (JTJ). In the special 
class of fields which only depend on one rectilinear coordinate n ■ x this equation can be 
Fourier transformed (three-dimensionally) into an ordinary differential equation: 



i(7 ■ n) 



d(n ■ x) 



+ 7 ■ k + 7 ■ A(n ■ x) — m 



Gh{ti • x,n ■ y,K) = 0, 



(2) 



with the conserved three-dimensional momentum coordinate k = k — nd(k ■ x)/d(n ■ x) 
orthogonal to n and where k stands for the four-momentum. As a further ansatz the matrix 
function Gn{n ■ x, n ■ y, k) is to be a functional of another matrix function gn{n ■ x,n ■ y, k) 
with a special property for the derivative: 



d 



d{n ■ x 



-G H [gH{n -x,n-y, k)] 



d(n ■ x 



-g H {n-x,n-y, k) 



G H [gH{n- x,n-y,K)], (3) 



which looks like the derivative of an exponential function, but is not quite due to the 
matrix structure. Provided a function gn(n • x,n ■ y, k) exists which satisfies equation Q, 
the form of the functional Gn[gH{n • x,n • y, k)] can be determined. Exploiting the above 
property leads to a factorisation in equation (J2J): 



£(7 ■ n) 



d 

d(n ■ x 



-9h[t l ' x i n ' Vi K ) + 7 ' K + 7 ' A(n ■ x) — in 



x G H [g H (n -x,n- y, «)] = 0. 



(4) 
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Given the existence of a solution Gh other than the trivial solution, its matrix structure 
can be inverted. Multiplication with the inverse of the solution from the right then yields a 
differential equation for the function gn'- 

z(7 • ^9n(n ■ x, n ■ y, k) + 7 • k + 7 • A(n ■ x) — m = 0. (5) 

This ordinary differential equation can be solved by direct integration where the initial 
condition <7#(n • x — n ■ y, n ■ y, k) — is chosen: 

7 • n f n ' x 

g H (n-x,n-y,K) = i — — / d(n ■ ^[7 ■ k + 7 • A(n ■ £) - m]. (6) 

Here it is necessary to require n 2 7^ 0. Otherwise, the matrix 7 • n does not posses an 
inverse because of det{7-n} = (n 2 ) 2 . For cases with n 2 = a different treatment is necessary. 

In general, the argument function gn does not commute with itself at different points 
n ■ x. This is not (only) due to non-Abelian charges which might be included in the vector 
field A but to the non-commutative nature of the elements of the Clifford algebra. Thus the 
solution of equation (jSJ) is not an exponential function but a path-ordered exponential: 

G H [gH(n -x,n-y, k)] = Vexp{g H (n -x,n-y, «)}. (7) 

A sufficient but not necessary condition for its existence is that the norm of the integrand 
in equation © is bounded. It has to be noted that the invariance of an integral under 
simultaneous exchange of the integration boundaries and inversion of the sign cannot be 
used in equation JHJ) because the path ordering would be reversed. 

If a more general initial condition had been chosen in equation (JOJ) the additional addend 
9H{n ■ y,n ■ y, k) would have lacked an ordering parameter necessary for the path-ordering. 
Hence, it could only be treated by always putting it to the rhs of the remaining path-ordered 
exponential. This would have led to an extra factor x exp{g#(n ■ y,n ■ y,k)}. As here a 
homogeneous differential equation is investigated this factor does not lead to independent 
solutions. 

Hereafter one has only to distinguish between the cases n 2 > and n 2 < 0, because every 
field A = A[n ■ x) can be transformed into a field A = A(n' ■ x) with sgn(n 2 ) = sgn(n /2 ) by a 
Lorentz transformation. Overall factors in front of the normal vectors can be absorbed in a 
redefinition of the vector potential A. So, for the sake of simplicity it suffices to investigate 
one special case per class of fields. This is going to be done for the cases of n M = (1,0,0,0) 
and n M = (0,0,0,1). In situations where n 2 = 0, rotations in three-space can turn any 
normal vector n into n M = (1, 0, 0, — 1)/\/2. 

A. Time-like coordinates 

For a purely time-dependent field, the solution for gH(xo,yo,k) in equation (jHJ) is given 

by: 

9h(x q , y , k) = 27 / d^ol^kj + 7 • A(£ ) - m], (8) 
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with j G {1,2,3}. Constructing the matrix function GH(x ,y ,k) by putting equation 
(JEJ) into equation (J2J) leads to: 

G H (x , yo, k) = V exp {27 jT° dftkfhj + 7 • A(£o) - m]} . (9) 

In the following, various approximations are studied in order to learn more about the 
above solution. 

1. Weak- field approximation 

It is useful to investigate the case of a vanishing gauge field A = 0. One then sees that 
the argument g°u{xQ — yo,k) now commutes with itself at different space-time points. The 
path-ordered exponential can now be replaced by an exponential function. The exponential 
function of matrices can be recast into exponential functions of scalar arguments multiplied 
with matrices: 

G H (*o-yo,k) = 7 ° 7 °" + ^ = m e+-^> + 7° 7 °" " P- + m e-"('°-">, (10) 



with to = y\k\ 2 + m 2 . 

Standard perturbation theory for small gauge fields A uo which can be interpreted as 
an ultraviolet approximation is obtained by expanding the exact solution in powers of A. 
Prior to this, it has to be rewritten in order to include all powers of the momenta and the 
mass with every factor of the field. The path-ordered exponential can be expressed as: 



N-l 




The interval [%,2/q] is decomposed into N disjoint pieces with the lengths A£g which 
need not be equal and each with an inner point t^f* . These are arranged according to 
x = Co < £1 < ••• < 6v < 6v+i = Vo for x < y or x = Co > 6 > ••• > 6v > 6v+i = Vo 
for xq > y$. V indicates that the factors are ordered with respect to the index v where the 
term with the lowest index is put furthest to the left. The expression can now be sorted 
with respect to powers of the field A: 

N N-l N-l N-l 

G H {x ,y ,k) = Jim^E EE-- E 

1=0 n 1= Ji2="-i+l ni=ni-i+l 

n\— 1 

II (l + H°^ L) Wk 3 -m]) x 

L=0 

U2 — 1 

x J] (l+*7°A^ L) [7^--m]) 

L=ni+1 



* 7 °7 " A(d ni) )A^ ni) ] 



*7°7 " A($»>)A& 
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X ... X 

n ; — 1 

L=n ; _ 1 + l 
N-l 

x II i 1 + i l°^o L) h i k j -m 

L=ni+1 



(12) 



For I = there are no further sums over rij. Sums and products are not taken into 
account if the starting index is greater than the ending index. In the limit iV — > oo the 
outer sum over the powers I of the gauge field A becomes an infinite sum, the intermediate 
sums turn into integrals over simplices, and the products give path-ordered exponentials. In 
fact, their arguments commute at every point, thus the path-ordering can be dropped here: 



G H (x ,y ,k) = 


oo 

1=0 J 


1 / d£ 2 ... / 

yo Jyo Jyo 






exp 


\il°[l 3 kj - m}(x - 


-6)}x 


X 


exp 


{z 7 V'%-™](£i- 


6)}x 


X 


... X 






X 


exp 




-6)} 


X 


exp 


{h^kj-m}^- 


yo)}- 



i 7 V^(£i) 
*7°7 " 

i 7 ° 7 ■ A(d) 



x 
x 



X 



(13) 



This expression is a uniformly and absolutely converging series representation for a path- 
ordered exponential: 

G H (x ,y , k) = G° H (xo ~ Vo, ^)Pexp || 10 d^ G° H (yo - Co, %7°7 • ^(£o)]G"^(£o - Jfo, £)} ■ 

(14) 

The above derivation is a special case of a more general identity for the type of path- 
ordered exponentials encountered here (see appendix A). The expansion of this formula in 
powers of A yields: 



— # p. — * 

G H (x ,y ,k) = G H (x -y ,k) + 

+ / d^G H (xo-^k)[t 1 ° 1 -A^ )}G° H ^o-yo,k) + 
Jyo 

r x o r£o -, -, 

+ diJ d Vo G° H (xo ~ fo, ^7°7 ■ A(^)]G° H ^o ~ Vo, k) x 
Jyo Jyo 

x[i 1 ° 1 -A{r ]0 )}G H {r} -y ,k) + 



+ ... 



(15) 



Up to now, only the solution GH(xo,yo, k) of the homogeneous Dirac equation in the 
mixed representation has been investigated. According to the equation of motion ((TJ the 



7 



inhomogeneous solution iG(x , y , k)^ must jump by one at x = y . The retarded propa- 
gator Gr(xo,Vo, k) vanishes for negative time differences xo — yo < 0. Due to the previous 
requirement on the argument gu(xo = yo,yo,k) = one has for the homogeneous solution 
Gh(x = yo,yo, k) = 1. Hence, the following condition has to be fulfilled in order to relate 
the latter and the retarded propagator: 

iG R (x , yo, £)7° = 9(x - yo)G H (x , y , k). (16) 

It should be noted that if the case of a different, more general coordinate with n 2 > 
should have been investigated at this point, the additional requirement n > would be 
needed here in oder to ensure that really the retarded propagator is obtained. However, this 
can always be achieved by a redefinition of the funcitonal form of the vector potential. 

All results obtained for the homogeneous solution of Dirac's equation in the present 
mixed representation are linked directly to the Green's function G R (x , y , k) by equation 
(Unj). After putting equation (JT3J) into the previous expression, the Heaviside function can 
be multiplied to every free homogeneous solution G° H (( , k): 

iG R (x ,y ,k)<y° = 9(x - y )G° H (x - y , k) + 

+ / d^eixo-^GUxo-^m^-A^o)} x 
Jyo 

^o-yo)G° H (^-yo,k) + 

+ f ° d£ [*° dri 9(x - £ )GP H (x - f , k) [i 7 ° 7 ■ A(£ )] x 
Jyo J yo 

x0(tio - ^o)G° H ^o - Vo, fc)[i7°7 • Mvo)] * 
*0(Vo ~ y )G° H {r]o - y , k) + 
+ ... . (17) 

This is possible due to the idempotency of the Heaviside function and the fact that 
9(xo — £o)#(£o - 2/o ) — @( x o — Vo) if £o G [xq, yo)- Subsequently, in accordance with equation 
([16)1 . the result can be reexpressed in terms of free Greens's functions: 

iG R (x ,y ,k) = iG R (x -y ,k) + 

+ / d^iG° R (x -^,k)[n-A^ Q )]iG R (^-y ,k) + 

+ / d£ / d Vo iG R (x - f , k)[i-y ■ A(£ )}iG R (£ - m , k) x 

x [«7 ■ A (Vo)]iG° R (r] - y , k) + 
+ ... . (18) 

Note that in the literature slightly different definitions exist for the propagator which 
account for the various occurrences of the imaginary unit i. The full retarded propagator 
Gr(xo, yo, k) inherits the asymptotic behaviour of the free retarded propagator G° R (xo—yo, k) 
by virtue of the above formula (fTSj) . 
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FIG. 1: Contour integration in the complex /co-plane for the determination of the Green's function 
with the correct assymptotic behaviour. The circles indicate the pairs of positions to which the 
poles are moved off the real axis by virtue of the corresponding e-prescription for the retarded 
(black) and the advanced (white) propagator. The squares show the position of the poles for the 
Feynman (white) and the reverse Feynman (black) propagators. For the retarded and the advanced 
propagators two poles or no pole is inside a given contour. For the Feynman and reverse Feynman 
propagators exactly one pole is always inside any contour. 



The full Feynman, i.e. time-ordered propagator cannot be expressed as a path-ordered 
exponential because it is defined with mixed boundary conditions: for the positive energy 
components at Xq — *> — oo and for the negative energy components at Xq — > +00. This can 
also be seen from the free Feynman propagator 



iG° F (x -yo,k) = e(x -y ) l0uJ ^ + m e -M*o-tt,) + 

+ % - Xo f° u + ^- m e ^o-yo) (19) 

which is a singular object in this and every mixed representation. That can be understood 
by looking at figure H Thus it is impossible to take its logarithm and express it as an 
exponential function. This is why the Feynman propagator cannot be equal to a path- 
ordered exponential of the form Q. 



Figure ^ shows the contour integrations in the complex /co-plane which have to be carried 
out in order to determine the contributions from the different poles of the corresponding 
propagator in momentum representation to that propagator in the mixed representation. 
Every pole included inside a contour results in an additive contribution to the free propagator 
proportional to one of the matrices 7 • k ± m. On shell, i.e. for k 2 = m 2 , these are singular. 
The circles in figure [T] belong to the retarded (black) and the advanced (white) propagators. 
It is important to note that either none of the poles is included in a contour or both. This 
means that, if one of these two propagators is non-zero, the two singular matrices occur in a 
non-trivial linear combination, which yields an invertible matrix structure. This is different 
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for the Feynman propagator and its relatives. There only one pole at a time is included in 
a contour. Thus these propagators are always non-invertible. 



2. Strong-field approximation 

The previous expansion which is appropriate for weak fields A(t) < w could be interpreted 
as an ultraviolet approximation. An infrared expansion requires a strong field A(t) ^> uj. 
It can be obtained by applying the resummation formula of appendix A in a different way. 
Resumming all scattering processes with the field for each power of the momentum one 
obtains: 



G H (x ,y ,k) = Pexp jijHdftTV^o)} x 



x Vexp 



H d£ Vexp \i [ V0 ^ o7 °7 • A(9 )\ {il°[l 3 k - m]} x 
Jyo < J to ) 

rto 



x Vexpli f ( ° ^ o7 °7 • A(6 
I Jyo 



(20) 



Now, one could start to expand the outer path-ordered exponential in powers of the 
momentum term: 



G H (x , y , k) = Vexp \i H ^o7°7 • Afo) \ + 

^ Jyo 

+ f d^VexpU I (lH ( f."~ ■ ,4(A>) } \h" y'k, - m\ \ x 
Jyo L J to 



d£ V exp |i £° ^o7°7 • M^o)] {i^l-fkj - m}} 
xVexp ji jf & ^o7°7 • A((9 )| + 

+ f ° ^0 / & d m Vexp U r ^o7°7 ' A(0 O )| {nV^ - ™]} x 
Jyo Jyo y J to > 

xV exp \i £° ^o7°7 • A((9 )| {i 7 V% - m]} x 

xPexp^£ ^ 7 7-A(^ )} + 
+ ... . *° (21) 

This corresponds to an expansion in powers of the on-shell energy u which can be un- 
derstood by noting that: {^°[^kj — m]) 2 = uj 2 . The weak-field approximation is based on 
the investigation of a given number of the - otherwise freely propagating - particle with 
the field which could be termed " accelerations" . The strong-field approach comes up to an 
expansion in powers of what could be called the " inertia" because uj equals the (asymptotic) 
relativistic mass. In the lowest order of the strong-field approximation the propagation of 
a particle without relativistic mass is governed only by the field. The higher order terms 
accord for deviations due to non- vanishing uj. 
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3. Abelian approximation 



All of the above approximations in form of an expansion with respect to some part of 
the exponent are based on equation (jA4j) . This is different for the Abelian approximation 
scheme(s), i.e. a commutative approximation with respect to the Clifford and the charge 
group algebra. The lowest order G^f(xo,yQ,k) of the Abelian approximation is given by 
omitting the path-ordering in equation (JHJ): 



Gi°(x , y , k) = exp j* 7 ° f° dtotfkj + 7 ■ A(£ ) - m]} . (22) 

Higher order approximations are not given by additive terms but by splitting equation 
dm at an ordered set of points: 



G^{x , y , k)=vf[ exp |z 7 ° t d^h + 7 • A{£ ) - rn}\ , (23) 

with zo = £o < £i < ... < £/v < 6v+i = Vo for x < y or x = £ > 6 > ••• > 6v > 
6v+i = Vo f° r > Vo- denotes that the factors are ordered with respect to the index v 
with the lowest index furthest to the left. The choice of the intermediate points is not 
unique, but in the limit of infinitely small intervals the result always becomes exact: 



lim GH N (x ,y ,k) 

N— >oo 

N 



f r& 

lim J] exp i 7 ° / d^kj + 7 ■ Afa) - m] 



(xo - yof 



N f 

= lim J] 1 + *7°(& - M b% + 7 ■ A(Co) - m] + 

N -*°°v^b { \ N 

= V exp («7° f ° d^T^i + 7 • A(fo) - m] 
= G H (x ,y ,k) (24) 

To estimate the error for an interval width of ?/o ~ = 2A, compare the lowest order 
result to the first order where the interval is divided into two halves exactly. 



AGi = G%° (0,2 A, k)-G% {0,2 A, k) = 

= Gjj°{0, 2 A, k) - Gj?{Q, A, k)G*°(A, 2 A, k) (25) 

With the help of the Baker-Campbell-Hausdorff formula: 



AG* = exp [g H (0, 2A, %)} - exp [g H (0, 2A, k) + [g H (0, A, k),g H (A, 2A, k)] + 0(A 4 )} = 

= -[g H (0,A,k),g H (A,2A,k)] + O(A A ) = 

= -[g H (0,A,k),g H (0,2A,k)] + O(A 4 ) = 

= -[g H (0,A,k),dg H (0,A,k)/dA]A + O(A 4 ). (26) 
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The first occurence of 0(A 4 ) results from a Taylor expansion of secondary and higher 
commutators. Thus, in leading order of the width of the interval A, the error is proportional 
to the commutator of the exponent gn and its first derivative at an intermediate point of 
the interval. For a constant integrand dgn(0,A,k)/dA, i.e. for a constant gauge field A, 
the Abelian approximation is exact. Higher order terms are required for fields that lead 
to a commutator [<?#(0, A, k), dgji(0, A, k)/dA]A not small against G H °(0, 2A, k). (This 
comparison must be based on the definition of an adequate norm.) 

This result can be compared to the error estimate for the standard form of express- 
ing a path-ordered exponential as a product of linear factors (see equation((TTJ)). For that 
representation one finds: 

AG H = 1 + 2g H (0, A, k)-[l + g H (0, A/2, k)} [1 + g H (0, 3A/2, k)} + 0(A 2 ) = 
= [2g H (0, A, k) - g H (0, A/2, k) - g H (0, 3A/2, k)} - 

-g H (0, A/2, k)g H {0, 3A/2, k) + 0{A 2 ) = 
= -g H (0, A, kf + [g H (0, A, k), dg H (0, A, k)/dA]A/2 + 0(A 2 ). (27) 

Contrary to AG H , AGh does not become zero for a constant gauge field. In leading order 
it depends on the actual value of the exponent gn- Thus, its convergence becomes slow not 
only for rapid changes of the gauge field but also for large values of the field and/or large 
energies. Even the free propagator then needs many terms to be approximated sufficiently 
well. 



B. Space-like coordinates 



The general solution scheme for a classical field depending on an arbitrary rectilinear 
coordinate n ■ x leading to equation (J2J) with the argument (jUJ) always yields a Green's 
function whose boundary conditions are given on a plane normal to n. Boundary conditions 
for propagators are given on surfaces with time-like normal vectors n 2 > 0. Hence, for a 
field only depending on the X3-coordinate {n^ = (0, 0, 0, 1)), only a Green's function but not 
a propagator is given by equations (JJJ) and (JHJ). 

This can also be seen directly. The solution for the argument gH(%3, V3', ko, kr) according 
to equation © is given by: 



gH{x3,y 3 ; k , k T ) = -*7 3 H d^°k + r J h + 7 ■ A(&) - m), (28) 

with an implicit sum over J G {1,2}. Repeating the steps that led to the free homoge- 
neous solution of the Dirac equation in the case n 2 > in equation (jTD|) yields: 



, i 3 -7 3 \/ fc o -m 2 T + 7% + J ] kj - m r 

G° H {x 3 -y 3 ;k ,k T ) = 7 3 * e V( 

2^J(k y-m: 



ko) 2 -m 2 T (x z -y:i) _|_ 



+ 7 d 



"7 3 \/ (^o) 2 — m T ~ 7°^o — l J kj + m 



2J{k Y-m 2 



(29) 
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with the transverse mass mr = y |&t| 2 + m 2 . This expression, multiplied with 8(x3 — y^j 
in order to obtain a Green's function from the homogeneous solution, is not proportional to 
the free retarded propagator in this mixed representation. 



C. Light-like coordinates 

As mentioned before, the present way to derive a homogeneous solution cannot be followed 
if the four-vector n is light-like, because in that case 7 • n has no inverse. However, for light- 
like coordinates there is a different approach that leads to a solution for Gh- In the case 
where the normal vector is = (1,0,0, -l)/y/2 equation © becomes: 



d 
dx- 



+ 7+[A;_ + - 7t • [hr + At(x~)\ + r y-A+(x-) - 



m 



G H (x-,y-; k-, k T ) 



0. 



(30) 



with v± = [vq ± v 3 ]/y/2 and v± = v T where v G {7, x, k, A(xJ)}. Noting that 7 + 7_/2 
and 7_7 + /2 are two projection operators which project into disjoint subspaces of the Clif- 
ford algebra and satisfy the completeness relation 7+7- + 7-7+ = 2 the matrix func- 
tion Gh(x-,V-] k_, k T ) can be split into 2Gh(x-, y~\ k_, k T ) = 7 + G_(x_, y_\ k_, hr) + 
7_G + (a;_, ?/_; hr) with G±(x~, y_; fey) = 7± &t). The argument will 

be suppressed in the following but until the end of the subsection the above mixed repre- 
sentation is addressed. Using this decomposition in equation (j3Uj) leads to: 

l ' t 2 ± ^ G - + 2j 2 h(A; - + A - )G+ ~ [1t ' ^ T + At) + m] ° H + ~ t 2 ±A+G - = °" (31) 

Use has been made of the idempotency of the projectors (7±7 T /2) 2 = 7±7 T /2 and their 
projection properties (7±7+/2)7 T = and (7±7+/2)7± = j±. From here, two equations can 
be obtained with the help of the projection operators: 



i ~d~x~°~ + A+G ~ ~~ \ ' ^ T + ^ + m ^- G+ = 
(k- + AJ)G+ - -[7r • {hr + A T ) + mfr+GL = 0. 



(32) 



The second equation is purely algebraic and can be used to replace G + in the first. 



i— G_ + - 

dx- 2 



7t ■ (kr + At) + m (k- + AJ) 7V ■ (hr + At) — m 



G- = 0. 



(33) 



When postulating a connection between the matrix function G_ and another g_ in direct 
analogy to equation (J3J) the resulting differential equation is given by: 



. d 1 

l dx~J- = "2 



7t • (kr + A T ) + ml (k- + AJ)~ l U T • (kr + A T 



m 



A, 



(34) 
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The equation can be solved by direct integration. As already argued before, in general the 
functional G- \gJ\ is given by the path-ordered exponential of its argument g_. In the present 
situation already the absence of non-Abelian charges turns it into an ordinary exponential 
because then it only contains the neutral element for multiplication of the Clifford algebra. 
The other component of the matrix function is given by the second of the equations (|32j). 
Finally a homogeneous solution of the differential equation (I3U|) has the form: 



G H (x-,y-; k_, k T ) 



7h 



:7-7+ 



[k- + A_0_)] 1 {7V • [k T + A T {x-)\ - m} 



x Pexp 



- {7V • [k T + I T (£_)] + m) + A-tf-)]- 1 {7V • [k T + A T (£-)} + A+(£_ 

'(35) 

If one tries to construct a propagator with the help of this homogeneous solution it can 
only be retarded or advanced in the light-like coordinate x_. Alternative approaches can be 
found in 0]. 

In the next chapter the production of fermion-antifermion pairs is described based on the 
results for the fermion propagator in a field that depends on a time-like coordinate. 



III. FERMION-ANTIFERMION PAIR PRODUCTION 



Here the results for the full propagator in an external field depending on one rectilinear 
time-like coordinate are applied to the problem of particle production due to vacuum po- 
larisation. First it is argued where such a propagator is of use in describing the physics of 
fermions in a heavy-ion collision. Second, a detailed comparison of the different approxima- 
tion schemes with the full solution for a given model field is presented. 

This calculation can be understood in a twofold way. On the one hand the field could be 
really an external field in the sense of the production of particles via vacuum polarisation. 
It is determined by the dynamics of the physical system without taking the back reaction of 
the particle creation into account. This field is used to calculate how many particles would 
be produced in its presence. This approach is justified if the process of particle production 
constitutes merely a small perturbation. Whether this condition is fulfilled has to be checked 
afterwards. 

On the other hand, the field could already be a self-consistent solution of a system of 
equations. For this solution for the classical field one would like to know how many particles 
were created in the process. In the present scenario such a set of equations would include 
the Yang-Mills equations with the expectation value for the current of produced fermions 
and antifermions (J u ) and an (initially present) external current J v ext : 



JLt + (J v ) 



(36) 
■ i \Afi, Ay\ 

in the adjoint representation. The expectation value for the current can be obtained from 
the causal propagator 22l |: 



with the gauge field and the corresponding field tensor T, 



d[j,Av duA^i ' 



14 



(r)~tr{YGc(x,x)}, 



(37) 



where the trace is only running over the matrices of the Clifford algebra and with the 
definition: Gc(x,x) = lim e ^ [Gc {x + ne,x) + Gc(x,x + ne)]/2 with n 2 > 0. Higher- 
order radiative corrections are suppressed by powers of the coupling constant which are not 
compensated by powers of the classical field. The causal propagator can be reexpressed as 
a linear combination of the retarded, the advanced, and the on-shell propagator: 

G c (x,y) = ±[G R (x,y) + G A (x,y) + G s (x,y)}. (38) 

— # 

The advanced propagator can be obtained from the homogeneous solution Gn(xo,yo, k) 
for the equation of motion for the Dirac Green's function by the relation: 



iG A (x , y , £)7° = -8(y - x )G H {x , y , k). (39) 

The on-shell propagator can be reexpressed in terms of the retarded and advanced one- 
particle scattering operators [Ti| : 

Gs(x,y) = J ^27r5(A; 2 -m 2 )G^(g)r fi (g,fc)(7-A; + m)r A (A;,p)GVp). (40) 

The scattering operators are defined according to equation below. In the present 
framework all required propagators are known as functionals of the classical gauge fields. 
Hence, equation (|3b[) with equation (|37jl constitutes an integro-differential equation for the 
classical gauge field. Its solution would yield the form of the field. The expectation value 
for the produced fermions and antifermions could be calculated from this field. 



Let us consider a model for the classical radiation-field in an ultra-relativistic heavy- 
ion collision. According to Bjorken jT^], the mid-rapidity region in a heavy-ion collision is 
characterised by boost-invariant quantities, i.e. boost-invariant along the beam direction. 
Let us consider a central collision in a symmetric system in the center-of-mass frame. For 
absolute values of the longitudinal coordinate \z\ smaller than the kinematic time t in this 
frame of reference the dependence on proper time r = y/t 2 — z 2 is approximated well by a 
dependence on the kinematic time t (see figure EJ). Most of the energy is deposited during 
t < t in close to the collision point at z = Hence, in good approximation, for t > ti n 
and \z\ < t in , an in general proper-time r dependent energy density can be reexpressed as 
an energy density depending on kinematic time t. Let the entire energy density be initially 
stored in an electric field of the form E v (t) (component along the hyperbolas). In the 
present approach, this fact is consistently approximated by a storage of the energy in E z (t). 
In temporal gauge (A = 0) or even in Lorentz gauge this is equivalent to a gauge field 
A z (t). For this form of gauge field, the retarded propagator has been derived in the previous 
section. 
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FIG. 2: Time dependence as approximation to the situation found in the central region of a boost- 
invariant system. The proper time r is constant on the hyperbolas. 

In order to proceed, one needs to know how to describe particle production based on a 
given propagator. With the Fourier transform of the retarded propagator: 

one can write the implicit definition of the corresponding one-particle scattering operator 
T in momentum space as: 

G R (q,p) = (2n)W\q - p)G° R (p) + G° R (q) x T(q,p) x G° R (p). (42) 
Explicitly, it is given by: 

T(q,p)=j-A(q-p) + I (^t ( ^« 7 • A{q - k)G R (k, 1)7 • A{1 - p). (43) 

In the following T always denotes the retarded one-particle scattering operator. The Born 
approximation, the expansion to lowest order in the fields, is given by: 

T B (q,p) = 1 .A{q-p). (44) 

This term is always contained in the scattering operator. The different approximation 
schemes discussed in the last section lead to differences in the remaining non-Bornian part 
in equation (j53J). In the presence of a purely time-dependent field the retarded one-particle 
scattering operator becomes: 

T(q,p) = {27T) 3 5 {3 \q-p) x 



x 



7 ■ A(q -p ) + / dxodyoe+^e-™^ ■ A(x )G R (x ,y ,p} 7 ■ A{y ) 



(27r) 3 6^(q-p)T( qo ,p ) = (2-K) 3 5®(q-$ \T B (q 0l p ) + T NB (q 0l p )] , (45) 
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with the non-Bornian part T NB (qo,po). Due to the occurence of the Dirac 5-distribution 
the conservation of the total three-momentum becomes obvious. In purely time-dependent 
fields the fermion-antifermion pairs are always produced in a back-to-back configuration. 
Terms of higher order in the gauge field - for instance in the weak-field expansion - can be 
obtained by replacing the full propagator in the scattering operator by terms from equation 
()18|) . Analogously, replacing the full propagator by various approximations leads to the 
corresponding approximations for the one-particle scattering operator. 

From the retarded one-particle scattering operator Tr the expectation value of produced 
pairs can be obtained: 



{H) = I 2(2^2(2^ \*® T *b> ~P>^ 2 ' ^ 

where a summation over the spin degrees of freedom of the unit spinors u(q) and v(p) is 

understood and where po = <J\p\ 2 + m 2 and go — \J\q\ 2 + m2 - F°r the special form of the 
scattering operator in spatially homogeneous purely time-dependent situations this simplifies 
to: 



(n) = / l«(Po, -p)T(po, -Po)v(po, +p) | 2 , (47) 

where use has been made of the relation [8^ 3 '(p — q)] 2 = VS^(p — q)/(2w) 3 . Carrying out 
the spin summation leads to: 



M^tm = tr ( T(po , _ po) /Po + 7^- ' V Tt _ po)7 . 7%-y R + m | 
V d 6 p { po po J 

In order to gain some insight into the behaviour of the differential expectation value (or 
momentum spectrum) and some information on the quality of the different approximations 
without having to solve the Yang-Mills equations ()36|) beforehand, the different formulae are 
going to be evaluated for a special choice of the field: 



A^t)=g^A in e~ t ' t °e(t). (49) 

Many other forms could have been taken. This choice was also inspired by a numerical 
study [19j which indicates that the field decays in a similar fashion. In any case, the actual 
form of the classical field has to be determined in a self-consistent calculation. 

For this field, the one-particle scattering operator in Born approximation (}4*4"|) is given 

by: 



T B (2u) = l3Amt ° . (50) 
v ; l + 2it uj v ; 

The following subchapters show the various approximations of the remaining part of the 
retarded one-particle scattering operator in the field P9*]l. 
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1. The weak- field approximation 



The lowest-order weak-field term of the full propagator is given by the free propagator 
(equation (fTBj) together with equation (fTOj)). Here, the non-Bornian part of the one-particle 
scattering operator is: 

' 2u 1+1 2uj ' ~ ' v 7 

with 

Ti /y = A? f°<feo [ X0 dyoe^^e^^-^e-^e-^ = -. ^t ) 2 /2 

Jo Jo [1 — icut \[l — i(uj ± LU)t \ 

(52) 

2. The strong-field approximation 

For any purely time-dependent field, the general expression for the homogeneous solution 
in the lowest-order strong-field approximation: 

G I H R (xo,yo,k)=Vexpi [i £° d£o7°7 • Afo) } (53) 

is not much simpler to evaluate than the exact solution. However, for a field of constant 
direction A^{t) = A^ x f(t) the path-ordering can be dropped. In the lowest-order strong- 
field approximation the non-Bornian part of the one-particle scattering operator is: 



T IR = _ ? 7^W T+ /« _ i t+^_ T IR^ (54) 



with: 



poo fZCO c 

rl R = A in 2 dxo dy e iuj{xo+Vo) e- xo/to e- yo/to exp[TiA in t [e- xo/to - e~ yo/ 

i a ., 2 ^ (±iAintal 1 ^ (± l A m t o r 1 f . 

- {A in t 2^ i n — ~~r 2^ i ; — rr, — <-,• . > \^ 

~ ill /i + 1 - ttot ^ i/! z/ + /i + 2- 2iut 

where use has been made of the uniform convergence of the exponential series for bounded 
arguments. With formula 6.5.29 in [2oj | : 



1 °° (— z ) n 

r(a,z) = —Y /7 ±- 2 V1 (56) 



for a bounded norm of Aj n t : 
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Tl R = (A m to) 2 E ^ Ai f o) " 1 , 7 *0u + 2 - 2z^ , ^Anto)r(/i + 2 - 2z^ ). 

(57) 

Formula 6.5.4 in |2C 



leads to: 



^ /i! 1 - 2Wt //i 

In the case of multiple charges, A in can be decomposed according to: A in = Af n T a where 
the T a are the generators of the corresponding algebra. Due to the requirement of unitarity 
these generators have to be hermitian. This is also true for any linear combination of the 
generators with real coefficients. Thus every matrix A^ n T a with real A^ n can be diagonalised, 
yielding: 

TV 

A a in T a = £ X n \n) (n\ , (60) 

n=l 

with the eigenvalues A n and the N orthonormal eigenvectors |n). The |n) (n| are projec- 
tors onto subspaces of different charges. Thus one gets: 

w = - e i«) h (TA„«„r<" f v^ +1 r 2 '?; T ' A " i ° ) - <«) 

n=l ix=\ V- 1 ~ ^0/^ 

5. The modified strong-field approximation 

In the present situation the special form of the field allows for a variation of the strong- 
field approximation, where the component of the momentum parallel to the field - in this 
case - is included in the exponent of the lowest-order expression A3 — » A3 + 

GW' (x ,y ,k) =Pexp jz £° d£o7°7 3 [*3 + A 3 (£ )]} ■ (62) 

In the lowest-order modified strong-field approximation the additional part of the one- 
particle scattering operator beyond the Born approximation is: 



T IR ' = -i^V - i£±£rl", (63) 
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with: 



Tl R = A in 2 / dx / dyoe^^^+^e-^e-^ exp [TiA in t [e- Xa ' tQ - e^* ]} . 

(64) 



Repeating the above steps leads to: 



A? 



1 7(/i + 1 - 2iut , =pA n t ) 



T± = ~ E \n) (n\ {T\nto? wto E V , —^ -r^, (g5) 

± h jkrf i-i(u±k 3 )t /ii 1 ; 

4- The Abelian approximation 

In the lowest-order Abelian approximation the interacting part of the retarded propagator 



is: 



T A = +iT A + iT A , (66) 

with 

T a = A J r dxQ r dyo ^(.o + ,o) e ^(xo-,o) e -.oAo e -,oAc 7°^ ± [7 J kj - 7 3 ^3 + m] 
Jo Jo 2fi 

(67) 



with the generalised energy = \jmT 2j rK 2 and the generalised momentum 
K 3 = k 3 + A in to(e- x ^ - e-»»/*o)/(xo - y )- 

Decomposition with respect to multiple charges leads to: 

AT 

T± = E l«> ("I V x 

n=l 

jo Jo 2f2 n 

(68) 

with the generalised energy Q n = \JrriT 2 + (K 3 ) n 2 and the generalised momentum 
(Ks) n = k 3 + A n to(e _:ro ^ — e _J, °/* () )/(a;o — yo) belonging to the respective eigenvalue A n . 
The modification to the longitudinal momentum is equal to the arithmetic average of the 
gauge field over the interval [x , y }. Hence the Abelian approximation can be interpreted as 
the description of the propagation of the fermions with their arithmetically averaged canon- 
ical momentum. In the weak-field expansion they are propagated with their asymptotic 
kinematic momentum. Higher orders in the Abelian approximation scheme make better 
approximations similar to a Fourier series. The particle is propagated with its canonical 
momentum averaged over every piece of the trajectory. The finer the partitioning of the 
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path, the closer the average canonical momentum is to its actual value in a particular inter- 
val. In the weak-field perturbative approximation scheme the particle is always propagated 
with its asymptotic kinetic momentum. For higher orders it only interacts with the field 
more and more often. 

For the strong-field approximation and the modified strong-field approximation the re- 
placement of the path-ordered exponential by an exponential function is only possible due to 
the special form of the field. In the Abelian approximation scheme this exchange is possible 
in the presence of an arbitrary field. 

The decomposition with respect to the charge projectors is also possible for the full 
non-Bornian part and the omnipresent Born part of the one-particle scattering operator 
T B = J2n=i \ n ) ( n \ %i ? - Hence, the whole operator can always be decomposed with respect 
to the same projectors T = E n =i \ n ) ( n \ \Xn + %f B }- m the squared expression needed to 
calculate the expectation value the contributions belonging to the different projectors do 
not mix. They lead to a sum over the expectation values for the different charge sub-spaces 
tr c {E£=i W) (n'|T„,£$ =1 |n"> (n" | T*„ } = ££=i |T„,| 2 tr c {\n') (n'\} = E^ =a |T„'| 2 , where 
tr c denotes the trace over the generators of the charge group. If the eigenvectors \n') are 
normalised, the remaining trace is equal to unity. Due to these facts, it suffices to compare 
the contributions from the different approximation schemes for one of the sub-spaces. 

It is always possible to measure all momenta, energies, and gauge field strengths in units 
of a scale parameter with the dimension of momentum. Then all lengths and times have 
to be given in units of inverse momenta. In the following, the eigenvalue belonging to the 
corresponding subspace is chosen as scale parameter and is going to be called A in /g again. 
The calculations are carried out assuming that all the energy of the system is included in 
one of the sub-spaces. 

The expected energy density produced in a central heavy-ion collision at LHC (Pb-Pb 
at y/s= 5.5 TeV) is e ~ 1000 GeV/fm 3 (§, 0, For the strong coupling constant one 
expects: a s ~ 0.15 [10]. If all the energy density was deposited in the field sector a rough 

estimate for the initial gauge field magnitude would be A in m \J g\f2e ~ 2GeV. For RHIC 
(Au-Au at \fs= 200 GeV) the typical coupling constant is around ot s ~ 0.33 and the initial 
energy density e ~ 50 GeV/fm 3 . This would lead to A in « lGeV. With decay times in the 
range from O.lfm/c to 0.5fm/c this leads to A in to between 0.5 and 5.0. Here only massless 
particles are investigated. 

The expressions for the Born (J50|) and the weak-field approximation (|51|) can be eval- 
uated straightforwardly. For the strong-field (fo^j) and the modified strong-field approach 
(f()3*|) the few first terms of the infinite series representations suffice for obtaining an accu- 
rate result. The integrals for the Abelian approximation have to be treated with standard 
numerical methods. The exact solution requires the handling of path-ordered exponentials 
and subsequent integrations. 



The general aspects of the exact solution for the momentum spectrum (J48|) are best 
seen in figures (J2K) and (@J). As a function of the transverse momentum it peaks once 
and shows no further relative extrema or other distinct structures. For increasing values 
of the parameter A in t (from figure (jl^) over figure (j2K) to figure (jib)) the peak in the 
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FIG. 3: Momentum spectrum of produced massless fermion-antifermion pairs versus transverse 
momentum and with the decay time toAi n = 2.0: a) Exact result for different values of the 
longitudinal momentum. In plots b), c), and d) the longitudinal momentum is fixed at = 0.1 Ai n 

b) Exact result (solid) compared to the Born (dashed) and the weak-field (grey) approximation. 

c) Exact result (solid) compared to the strong-field (dashed) and the modified strong-field (grey) 
approximation, d) Exact result (solid) compared to the Abelian (dashed) approximation. 



transverse momentum spectrum becomes more pronounced for a fixed value of the longi- 
tudinal momentum k^. In other words it increases in height and decreases in width (see 
especially the different scale of the transverse momentum axis in figure (jib)). Actually, 
the differential expectation value is a function of the variable ooto. Hence the width of 
the transverse-momentum spectrum for massless particles at mid-rapidity scales exactly in- 
versely proportionally to t . The same holds still after uj and to have been rescaled with 
A in . For fields of a functional form analogous to that of the present special model field the 
peak height seems to be strictly monotonically decreasing with increasing longitudinal mo- 
mentum, as is suggested in figure (J3k). Further, at zero transverse momentum no particles 
are produced; the fermions and antifermions are never produced with momenta along the 
direction of the field but preferentially with momenta perpendicular to the field. 

A comparison of the different approaches shows that for large momenta all approximations 
and the exact solution tend towards the Born result. This is due to the form of the one- 
particle scattering operator (j%5|) . Together with the Born approximation all other graphs 
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FIG. 4: Momentum spectrum (solid) of produced massless fermion-antifermion pairs versus trans- 
verse momentum compared to the Born (dashed) and the Abelian (grey) approximation for fixed 
longitudinal momentum k^, = 0.1 Ai n and for different values of the decay time: a) Ai n to = 0.5 and 
b) Ai n to = 4.0. In plot b) the enhanced strong-field approximation (dashed grey) is shown, too. 

tend toward zero for higher particle energies. As shown in figure (JHb) the Born approach 
overestimates the exact value for low momenta but underestimates it for high momenta. The 
weak-field approximation is an improvement compared to the Born approach for most values 
of the transverse momentum. Looking at figure (J2b) the strong-field and the modified strong- 
field are generally closer to the exact result than the weak-field approximation. However, 
for more general forms of time-dependent fields the propagators in these schemes are not 
much simpler to deal with than the full one. The modified strong-field approximation even 
ceases to be available because the terms longitudinal and transverse might no longer be 
well-defined with respect to the field. For all momenta the Abelian approximation scheme 
(see figure (jUli)) is closest to the exact values. The largest deviations are found for small 
energies and large values of the parameter A in t (compare figures (jUli), (jlK)) and (El 3 ))- The 
reason is that there the situation is maximally non-Abelian, i.e. there the condition that 
the typical commutator of the exponent gni^o, yo, k) at different points is negligible with 
respect to the typical propagator is least well satisfied (see also equation (J2rJ|) ). While for 
low values of A in to the Born approximation is reasonably good it is not appropriate for large 
values (see figure (JU)). 

IV. SUMMARY 

The exact homogeneous solutions for the Dirac equation in a gauge field depending on 
one rectilinear coordinate has been presented. An alternative way had to be taken for a 
dependence on a light-like coordinate. In the case where this coordinate was time-like, the 
retarded propagator has been constructed from the homogeneous solution. The analogous 
result for a space-like coordinate was seen to constitute a Dirac Green's function but not a 
propagator. 

For the situation of a time-like coordinate various approximation schemes for the exact 
solution have been determined. Explicitly, these are the weak-field approximation, the 
strong-field approach, and the Abelian approximation. Additionally, a larger variety of 
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approximations can be obtained with the help of the general resummation formula. 

Subsequently, the retarded fermion propagator and all the lowest orders of the various 
approximation schemes in the presence of a gauge field depending on one rectilinear time- 
like coordinate have been used to calculate the momentum spectrum of produced fermion- 
antifermion pairs. The resulting expressions are evaluated for a decaying model field and 
the results are mutually compared for parameters expected to be found in ultra-relativistic 
heavy-ion collisions. In this situation an additional modified strong-field approximation 
could be obtained. 

In the present situation, the exact momentum spectrum is a singly peaked function of the 
transverse momentum with no further distinct structure. The quality of the approximations 
increases from the Born approach over the lowest-order weak-field, strong-field, modified 
strong-field, towards the Abelian approximation. It should be mentioned that in more 
general situations the strong-field and the modified strong-field approach are not much 
simpler to evaluate than the full result. The model parameter is A in t Q . It is the product of 
the initial magnitude of the gauge field A in and the decay time scale of the field. For the 
smallest expected values the Born approximation is still acceptable. Nevertheless, the other 
schemes like the Abelian or the enhanced strong-field are even better. For the highest values 
of the decay time only the latter come close to the exact result. Hence, in order to ensure the 
maximum possible independence from the scale parameter Ai n to without having to evaluate 
the exact solution it would be best to use the Abelian approximation for self-consistent 
calculations. 
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APPENDIX A: GENERAL RESUMMATION FORMULA 

In general, a path-ordered exponential with an integrand depending on a single variable 
£ can be rewritten in the following way: 

Vexpl df [5(fo)+C(fo)] = £ / d£ 2 ... 
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Making use of the group property valid for the present path-ordered exponentials: 

Pexp |£° dCoB(Co)) x Pexp |£° d£ £(£o)} = P exp |£° ^ 5(£o)} (A2) 
the above equation can be reexpressed as: 
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(A3) 



Again, by virtue of the group property the most compact form is given by: 
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Summarising, more general resummations are possible in which the dominant quantity 
can be chosen arbitrarily. Further, the above steps can be repeated so as to resum the 
obtained quantity several times, e.g. after splitting C(£o) into a dominant part and a 
deviation. 

The above result can also be obtained in another way. Write the Fourier transformed 

— * 

Dirac equation in a generic form, where the dependence on k will not be denoted in the 
following: 



d 
dxo 



B(x ) - C(x ) 



G H (x ,yo) = 0. 



With the ansatz 



(A5) 



G H (x ,yo) = U(x ,yo)G H (xo,yo) 
and the product rule for differentiation one obtains: 



(A6) 



—G H (x ,y ) 
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- [B(x ) + C(x )] U (x , y )G H (x , y Q ) = 0. 
When postulating that U(xo,yo) satisfies the differential equation 
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the above expression reduces to: 
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This formula can be transformed by multiplying with U 1 (xq,?/o) from the left: 



- U \x , y )C{x )U{x , y )G H (x , y ) = 0. (A10) 
As has been shown above, the solutions for the last equation and equation (jA8|) are given 

by 

6{x , y ) = V exp { £° d^U~\^ y )C(Z )U(£ , y )] (All) 

and 

U(x ,y ) =Vexp{£dZoB(£ )}, (A12) 

respectively. Making use of the group property (jA2|) the last line of equation (jA3|) is 
reproduced. These derivations can be repeated for any variable n ■ x. Only for n 2 = the 
Dirac equation cannot be written in the generic form (JA5j) . 
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